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Abstract 

Motivated by the recent experimental data, we have revisited the B — > ttK, tttt decays 
in the framework of QCD factorization, with inclusion of the important strong penguin 
corrections of order induced hy b ^ Dg*g* (D = d or s and g* denotes an off-shell 
gluon) transitions. We find that these higher order strong penguin contributions can 
provide ~ 30% enhancement to the penguin-dominated B ttK decay rates, and such 
an enhancement can improve the consistency between the theoretical predictions and the 
experimental data significantly, while for the tree-dominated B — > vrvr decays, these higher 
order contributions play only a minor role. When these strong penguin contributions are 
summed, only a small strong phase remains and the direct CP asymmetries get small 
corrections. We also find patterns of the ratios between the CP-averaged branching 
fractions remain nearly unaffected even after including these higher order corrections and 
the ttK puzzle still persists. Our results may indicate that to resolve the puzzle one would 
have to resort to new physics contributions in the electroweak penguin sector as found by 
Buras et al. 
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1 Introduction 



The study of exclusive hadronic 5-meson decays can provide not only an interesting avenue to 
understand the CP violation and flavor mixing of the quark sector in the Standard Model (SM), 
but also powerful means to probe different new physics scenarios beyond the SM. With the 
operation of i?-factory experiments, large amount of experimental data on hadronic 5-meson 
decays are being collected and measurements of previously known observables are becoming 
more and more precise. Thus, studies of the hadronic 5-meson decays have entered a precision 
era. 

With respect to the theoretical aspect, several novel methods have also been proposed to 
study exclusive hadronic B decays, such as the "naive" factorization (NF) [1] , the perturbative 
QCD method (pQCD) [2], the QCD factorization (QCDF) [3, 4], the soft collinear effective 
theory (SCET) [5], and so on. For quite a long time, the decay amplitudes for exclusive two- 
body hadronic B decays were estimated in the NF approach, and in many cases, this approach 
could provide the correct order of the magnitude of the branching fractions. However, it cannot 
predict the direct CP asymmetries properly due to the assumption of no strong rescattering in 
the final states. It is therefore no longer adequate to account for the new i?-factory data. The 
other methods mentioned above are proposed to supersede this conventional approach. Since 
we shall use QCDF approach in this paper, we would only focus on this approach below. 

The essence of the QCDF approach can be summarized as follows: since the h quark mass is 
much larger than the strong interaction scale Aqcd, in the heavy quark limit nib ^ ^qcd, the 
hadronic matrix elements relevant to two-body hadronic 5-meson decays can be represented 
in the factorization form [3] 



where Qi is the local four-quark operator in the effective weak Hamiltonian, ji 2 are bilinear 
quark currents, and Mi is the meson that picks up the spectator quark from the B meson, 
while M2 is the one that can be factored out from the {B,Mi) system. This scheme has 
incorporated elements of the NF approach (as the leading contribution) and the hard-scattering 
approach (as the sub- leading corrections). It provides a means to compute the hadronic matrix 
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elements systematically. In particular, the final-state strong interaction phases, which are very 
important for studying CP violation in 5-meson decays, are calculable from first principles with 
this formalism. Its accuracy is limited only by higher order power corrections to the heavy- 
quark limit and the uncertainties of theoretical input parameters such as quark masses, form 
factors, and the light-cone distribution amplitudes. Details about the conceptual foundations 
and the arguments of this approach could be found in Ref. [3, 4]. 

Among the two-body hadronic i?-meson decays, the charmless B ttK and B ^ titt 
modes are very interesting, since a significant interference of tree and penguin amplitudes is 
expected, and hence have been studied most extensively. Experimentally, all the four decay 
channels for B ^ nK {B^ tt^K^, B^ n^K^, B^ n^K^, and 5° ^ K^tt^) and the 
three ones for B ^ ttti {B^ ti^tx^ ^ B^ tt'^tt", and B^ 7r°7r°) have been observed with 
the CP-averaged branching ratios measured within a few percent errors by the CLEO [6, 7, 8], 
BaBar [9], and Belle [10] collaborations. The CP asymmetries in these decay modes have 
also been measured recently [11, 12, 13, 14, 15, 16, 17, 18, 19]. In particular, measurements 
of the direct CP asymmetry in B^ n^K^ have been recently achieved at the 5.7(j level by 
BaBar [13, 14] and Belle [15, 16, 17, 20]. All these experimental data can therefore provide very 
useful information for improving the existing model calculations. On the theoretical side, these 
decay modes have also been analyzed in detail within the QCDF formalism [21, 22, 23, 24]. Due 
to lack of precise experimental data at that time, no large discrepancies between the theoretical 
predictions and the experimental data were found. However, the current new 5-factory data for 
B ttK, tttt decays indicate some potential inconsistencies with the predictions based on this 
scheme. For example, new experimental data for B^ n^K^, vr^vr^ decay rates are significantly 
larger than the theoretical predictions with this scheme. In addition, predictions for the direct 
CP asymmetries in these modes are also inconsistent with the data, even with the opposite sign 
for some processes [21, 25]. Moreover, the experimental results of the following ratios between 
the CP-averaged branching fractions for B ivK, vrvr decays [26, 27] 
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with numerical results compiled by the Heavy Flavor Averaging Group (HFAG) [28] , have shown 
very puzzling patterns [29, 30]. Within the SM, predictions based on the QCDF approach give 
Rc ~ Rni while the value for R is quite consistent with the experimental data [25]. The 
central values for and -Rqo calculated with the QCD factorization [25] give = 1.24 
and -Rqo = 0.07 as emphasized by Buras et al. [30], which are also inconsistent with the current 
experimental data. Though none of these exciting results is conclusive at the moment due 
to large uncertainties both theoretically and experimentally, it is important and interesting to 
take them seriously and to find out possible origins of these discrepancies. Recently, quite a 
lot of works have been done to study the implications of these new experimental data [30, 31, 
32, 33, 34, 35, 36, 37, 38, 39, 40, 41]. In this paper, we restrict ourselves to the possibility 
that these deviations result from our insufficient understanding of the hadronic dynamics and 
investigate the higher order strong penguin effects induced by 6 — > Dg*g* transitions, where 
D = d OT s, depends on the specific decay modes. The off-shell gluons g* are either emitted 
from the internal quark loops, external quark lines, or splitted off the virtual gluon of the strong 
penguin. 

As shown in literature [42, 43, 44, 45, 46, 47], contributions of the higher order b sgg pro- 
cess to the inclusive and semi-inclusive decay rates of 5-meson decays could be large compared 
to b ^ sg process. For example, in [45] , Greub and Liniger have found that the next-to-leading 
logarithmic result of B^^^{b sg) = (5.0 ±1.0) x 10"^ is more than a factor of two larger 
than the leading logarithmic one B^^{b sg) = (2.2 ± 0.8) x 10~^. In addition, in [47], we 
have found the higher order strong penguin could give large corrections to i? ^ (pXg. We also 
note that the large higher order chromo-magnetic penguin contributions have also been found 
by Mishima and Sanda[48] in PQCD factorization framework. Since the B tiK decays are 
dominated by strong penguin contributions, it is interesting to investigate these higher order 
b sg*g* strong penguin effects on these penguin-dominated processes. However, for self- 
consistent, we will also investigate these effects on the tree-dominated B ^ tttt decays. After 
direct calculations, we find that these higher order strong penguin contributions can provide 
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~ 30% enhancement to the penguin-dominated B — > vr/i decay rates, and such an enhancement 
can improve the consistency between the theoretical predictions and the experimental data ef- 
fectively. For tree-dominated B ^ im decays, however, their effects are quite small. Since 
the b Dg*g* strong penguin contributions contain only a relatively small strong phase, their 
effects on the direct CP asymmetries are also small. In addition, the patterns of the quantities 
R, Rc, Rni R+-1 and i?oo defined above remain unaffected even with these new contributions 
included. 

This paper is organized as follows: In Sec. 2, using the QCDF approach, we first calculate 
the B 71 K, TTTT decay amplitudes at the next-to- leading order in a^, and then take into 
account the b Dg*g* strong penguin contributions to the decay amplitudes. In Sec. 3, after 
presenting the theoretical input parameters relevant to our analysis, we give our numerical 
results for B txK and B ^ tttt decays. Some discussions on these higher order corrections 
and the 7 dependence of the relevant quantities are also presented. Finally, we conclude with a 
summary in Sec. 4. In Appendix A, we present the correction functions at next-to-leading order 
in Os- Explicit form for the quark loop functions are given in Appendix B. 



2 Decay amplitudes for B ttK^ tttt decays in QCDF 
approach 

2.1 The effective weak Hamiltonian for hadronic B decays 

In phenomenological treatment of the hadronic 5-meson decays, the starting point is the ef- 
fective weak Hamiltonian at the low energy [49, 50], which is obtained by integrating out the 
heavy degree of freedom (e.g. the top quark, and Z bosons in the SM) from the Lagrangian 
of the full theory. After using the unitarity relation — Aj = + Ac, it can be written as 

^eff = % E >^P (Cl + C2 + E Q^ + Q77 + Qsg) + h.C. , (7) 

V 2 p=u,c \ i=3,...,10 ^ 

where Ap = Vpb V*g (for b ^ s transition) and A^ = Vpb V*^ (for b ^ d transition) are products of 
the CKM matrix elements. The effective operators Qi govern a given decay process and their 
explicit form can be read as follows. 
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Current-current operators: 



{ph)v-A{Dp)v-A, 



Ql = {Vibj)v-A{DjPi)v-A 



- QCD-penguin operators: 

Qs = (£'6)y_A^g(gg)y-A, 
Qs = {Dh)v-A^q{.qq)v+A, 

- Electroweak penguin operators: 

Qi = {Db)v-AY.g¥i(^(l)v+A^ 
Qg = {Dh)y^A^qleq{qq)v-A, 

- Electro- and chromo-magnetic dipole operators: 



Qi = {Dibj)v-A {(lj(li)v-A , 
Qq = {Dibj)v-A ^q {qj<li)v+A , 
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Qs = {Dibj)v-A^q |eq(gjgj)y+A , 
Qio = {Dibj)v^A J2<i ^'^qiQj1i)v-A , (10) 



where {qiq2)v±A = 15)12, hj are colour indices, Cq are the electric charges of the quarks 

in units of |e|, and a summation over q = u,d, s,c,b is implied. For decay modes induced by 
the quark level b ^ d transition, D = d, while for 6 — s transition, D = s. 

The Wilson coefficients Cj(/i) in Eq. (7) represent all the contributions from physics with 
scale higher than /i ~ 0{mb) and have been reliably evaluated up to the next-to-leading loga- 
rithmic order. Numerical results for these coefficients evaluated at different scales can be found 
in [49]. 



2.2 Decay amplitudes at the next-to-leading order in ag 

Using the weak effective Hamiltonian given by Eq (7), we can now write the decay amplitudes 
for the general two-body hadronic B M1M2 decays as 
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(12) 



Then, the most essential theoretical problem obstructing the calculation of the hadronic B- 
meson decay amplitudes resides in the evaluation of the hadronic matrix elements of the local 
operators {MiM2\Q^\B) . Within the formalism of the QCDF, this quantity could be simplified 



Figure 1: Graphical representation of the factorization formula. Only one of the two 
form- factor terms in (13) is shown for simplicity. 

greatly in the heavy-quark limit. To leading power in Aqcd / "^b, but to all orders in perturbation 
theory, it obeys the following factorization formula [21] 

+ Tl^ * $B * $A/i * $Af2 , (13) 

where $m is the leading-twist light-cone distribution amplitude of the meson M, and the * 
products indicate an integration over the light-cone momentum fractions of the constituent 
quarks inside the mesons. The quantity F^^^^ denotes the B ^ M transition form factor. 
This formula is illustrated by the graphs shown in Fig. 1. 

In Eq. (13), the hard-scattering kernels T\,j^j and T}^ are calculable order by order with the 
perturbation theory. T\i^j starts at tree level, and at higher order in contains the "non- 
factorizable" corrections from the hard gluon exchange and the light-quark loops (penguin 
topologies). The hard "non-factorizable" interactions involving the spectator quark are part 
of the kernel T}^. At the leading order, T\.i ij = 1, T}^ = 0, and the QCDF formula reproduce 
the NF results. Nonperturbative effects are either suppressed by h^Qco/^b or parameterized 
in terms of the meson decay constants, the transition form factors F^^^'\ and the light-cone 
distribution amplitudes ^m- The relevant Feynman diagrams contributing to these kernels 
at the next-to- leading in are shown in Fig. 2. 

According to the arguments in [3], the weak annihilation contributions to the decay ampli- 
tudes are power suppressed compared to the leading spectator interaction in the heavy quark 
limit, and hence do not appear in the factorization formula (13). Nevertheless, as emphasized 
in [2, 51, 52], these contributions may be numerically important for realistic i?- meson decays. 
In particular, the annihilation contributions with QCD corrections could give potentially large 
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Figure 2: Order corrections to the hard-scattering kernels Tj^j - (coming from the 
diagrams (a)-(f)) and T/^ (from the last two diagrams). 




(a) (6) '"i (c) (d) 

Figure 3: The annihilation diagrams of order as- 

strong phases, hence large CP violation could be expected [2, 51]. It is therefore necessary 
to take these annihilation contributions into account. At leading order in a^, the annihilation 
kernels arise from the four diagrams shown in Fig. 3. They result in a further contribution to 
the hard-scattering kernel T/^ in the factorization formula. 

As indicated in the factorization formula (13), the meson light-cone distribution ampli- 
tudes (LCDAs) play an important role in the QCDF formalism. For convenience, we list the 
relevant formula as follows (details can be found in [53]) 

- LCDAs for B meson. In the heavy quark limit, the light-cone projector for the B meson 
in the momentum space can be expressed as [3, 53, 54] 

ifsfTiB 



with the normalization condition 



:i+?^)75 K(o+ 
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where ^ is the momentum fraction of the spectator quark in the B meson. For simplicity, we 
consider only the leading twist $f (0 contribution in this paper. Since almost all the momentum 
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of the B meson is carried by the heavy b quark, we expect tliat $f (0 = C>{^b/ ^qcd) and 

e = 0(AQCD/m,). 

- LCDAs for light mesons. For the hght-cone projector of hght pseudoscalar mesons in 
momentum space, we use the form given by [55] 
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where fp and p are the decay constant and the momentum of the meson. The parameter /ip = 
rn?p / imi^jj) +m2(/i)), with mi,2(At) being the current quark mass of the meson constituents, is 
proportional to the chiral quark condensate. $(x) is the leading- twist distribution amplitude, 
whereas ^p{x) the sub-leading twist (twist-3) one. All of them are normalized to 1. The quark 
and anti-quark momenta of meson constituents, ki and /c2; are defined respectively by 

P P 
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2xp ■ p 2[1 — x)p ■ p 

where p is a light-like vector whose 3-components point into the opposite direction of p. It is 
understood that only after the factor ki ■ k2 in the denominator of Eq. (16) cancelled, can we 
take the coUinear approximation, i.e., the momentum ki and /c2 can be set to xp and (1 — x) p, 
respectively. 

From now on, we denote by u the longitudinal momentum fraction of the constituent quark 
in the emitted meson M2, which can be factored out from the (5, Mi) system, and by v the 
momentum fraction of the quark in the recoiled meson Mi, which picks up the spectator quark 
from the decaying B meson. For B meson decaying into two light energetic hadronic final 
states, we define the light-cone distribution amplitudes by choosing the -|- direction along the 
decay path of the emission meson M2. 

Equipped with these necessary preliminaries, the four B ttK and the three B 
decay amplitudes can be expressed as [21, 25] 
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where the "chirally-enhanced" factor = r^^{fi) associated with the coefficients ag and ag is 



defined by 



r-(/i) = — 



(20) 



mb(/i) {mu,d{lJ,) + m^(/i)) ' 
with friqlfi) being the current quark mass and depending on the scale /i. The CP-conjugated 
decay amphtudes are obtained from the above expressions by just replacing A^'^ with A^'-**. 

In Eq. (18) and (19), we have defined x^'^^^^'^^^^ as the factorized amplitude with the meson 
M2 being factored out from the (5, Mi) system 



(21) 



In term of the decay constant and the transition form factors defined by [53, 56] 
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the factorized amplitude can be written as 

72 



(24) 



where we have combined the factor ^ in the effective Hamiltonian. The quantity x^^''^^^'^^^ 
associated with the annihilation coefficient bi and fc^"' is given by 

y(B,MiM2) _ a^F f f f 

— ^ JB J Ah JM2 



(25) 



The parameters = aj(MiM2) in Eq. (18) and (19) encode all the "non-factorizable" 
corrections up to next-to- leading order in a^, and are calculable with perturbative theory. The 
general form of these coefficients af can be written as [25] 
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where Cp = {N^ — l)/(2A''c), and A'",:; = 3 is the number of colors. The upper (lower) signs 
apply when i is odd (even) and the superscript 'p' should be omitted for i = 1,2. The first part 
in Eq. (26) corresponds to the NF results, and the remaining ones to the corrections up to the 
next-to- leading order in a^. The quantities Vi{M2) account for the one- loop vertex corrections, 
Hi^MiM^) for the hard spectator interactions, and Pf (M1M2) for the penguin contractions. In 
general, these quantities can be written as the convolution of the hard-scattering kernels with 
the meson distribution amplitudes. Explicit form for these quantities are relegated to Appendix 
A. 

The parameters bi = bi{MiM2) in Eq. (18) and (19) correspond to the weak annihilation 
contributions and are given as [25] 
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(27) 
(28) 
(29) 
(30) 



where we have omitted the argument "M1M2". These coefficients correspond to the current- 
current annihilation (61, ^2); the penguin annihilation (63,64), and the electro-weak penguin 
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annihilation (^3"', b^^), respectively. The explicit form for the building blocks A]^-' can be found 
in Appendix A. 

It should be noted that within the QCDF framework, all the nonfactorizable power sup- 
pressed contributions except for the hard spectator and the annihilation contributions are 
neglected. We have re-derived the above next-to-leading order formulas calculated by Beneke 
and Neubert [25], for which no deviation has been found. 

2.3 The b Dg*g* strong penguin contributions to the B ttKjTttt 
decays 

From the previous subsection, we can see that, up to next-to-leading order in as and to leading 
power in Aqcd/?^;), the strong-interaction phases originate from the imaginary parts of the 
functions g{u) and G{s,u), as defined in Eq. (70) and Eq. (76), respectively. The presence of 
strong-interaction phase in the penguin function G{s,u) is well known and commonly referred 
to the Bander-Silverman-Soni (ESS) mechanism [57]. The reliable calculation of the imaginary 
part of function g{u) arising from the hard gluon exchanging between the two outgoing mesons 
is a new product of the QCDF approach. However, recent experimental data indicate that 
there may exist extra new strong interaction phases in hadronic 5-meson decays. Since the 
b sgg transitions play an important role in the inclusive and semi-inclusive 5-meson decays 
as discussed in literature [42, 43, 44, 45, 46, 47, 48], in this section we shall generalize these 
results to exclusive two-body hadronic B decays, and investigate these b Dg*g* strong 
penguin contributions io B ^ ttK, vrvr decays. 

At the quark level, the b — > Dg*g* transitions can occur in many different manners as de- 
picted by Figs. 4-6. For example, one of the gluons can radiate from the external quark line, 
while the other one coming from the chro mo-magnetic dipole operator O^g as in Figs. 5(b) 
and 5(c) or from the internal quark loop in the QCD penguin diagrams in Figs. 6(b) and 6(c). 
On the other hand, the two gluons can also radiate from the internal quark loops in Figs. 6(d) 
and 6(e) or split off the virtual gluon of the strong penguin processes as shown by Figs. 5(a) 
and 6(a). Here we do not consider the diagrams of the category in Fig. 4, since their contri- 
butions can be absorbed into the definition of the B Mi transition form factors Figs. 4(a) 

2 

and 4(b) or further suppressed by y|^- It is easy to clarify this point by comparing the strengths 

12 



of Fig. 4(c) to that of Fig. 5(a). 



M2 M2 M2 




(d) (e) 



Figure 4: Representative diagrams induced by 6 ^ Dg*g* transition which are not evaluated. Here 
we give only the chromo-magnetic dipole operator Qgg contributions. With Ogg replaced by the other 
operators, the corresponding diagrams for these operators can also be obtained. 
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Figure 5: Chromo-magnetic operator Q^g contributions induced by 6 ^ Dg*g* transition. 




13 



As can be seen from Figs. 5 and 6, these penguin diagrams should be the dominant con- 
tributions of order a^, since they are not two-loop QCD diagrams and there is no additional 
suppression factor compared to the genuine two-loop contributions of order a^. Studies 
of these contributions could be helpful for understanding the higher order perturbative cor- 
rections within the QCDF formalism. In followings, we first discuss these higher order strong 
penguin contributions to decay modes with two light pseudoscalar mesons in the final states, 
B —>■ M1M2, and then specialize this general case to the B vrK, tttt decays and investigate 
the effect of these higher order corrections on the branching ratios and CP asymmetries for 
these modes. 

We start with the calculation of the diagrams in Fig. 5. In this case, the weak decay is 
induced by the chromo-magnetic dipole operator Osg- The calculation is straightforward with 
the result given by 
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(31) 



where At = Vf^V^* (for b ^ s transition) and A'^ = VtbV^*^ (for b ^ d transition) are products 
of the CKM matrix elements. As always, $a/ and denote the leading-twist and twist-3 
LCDAs of the pseudoscalar meson M in the final state, respectively. 

In calculation of the Feynman diagrams of Fig. 6, we follow the method proposed by Greub 
and Liniger [45]. First, we calculate the fermion loops in these individual diagrams, and then 
insert these building blocks into the entire diagrams to obtain the total contributions. In 
evaluating the internal quark loop diagrams, we shall adopt the naive dimensional regulariza- 
tion (NDR) scheme and the modified minimal subtraction (MS) scheme. In addition, we shall 
adopt the ad hoc Feynman gauge throughout this paper. Similar to the calculation of the 
penguin contractions in Appendix A, we should consider the two distinct contractions in the 
weak interaction vertex of these penguin diagrams. 



14 



Ql.3 Q4.6 




Figure 7: Building blocks I^{k) (associated with the contraction of the operators Qi^s) and I^{k) 
(associated with the contraction of the operators Q^^q) for Figs. 6(a)-6(c). 



As can be seen from Fig. 6, the first three diagrams have the same building block /^(/c) 
(corresponding to the contraction of operators Qi^s) or I^{k) (associated with the contractions 
of the operators Qi^). These building blocks are shown in Fig. 7 and given by 
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where k and T" is the momentum and the color generator of the off-shell gluon, Qs is the strong 
coupling constant, and the pole mass of the quark propagating in the quark loops. The free 
indices ^ and a should be contracted with the gluon propagator when inserting these building 
blocks into the entire diagrams. Here we have used the d dimension space-time as 0? = 4 — e. 
After performing the subtraction with the MS scheme, we get 
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(34) 
(35) 



with the function G{s,u) defined by Eq. (76). 

The sum of the fermion loops in the last two diagrams in Fig. 6 are denoted by the building 
block J^l{k,p) (corresponding to the contraction of operators Qi.a) or J^l{k,p) (corresponding 
to the contraction of operators Q^^), as depicted by Fig. 8. Using the decomposition advocated 
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by [44, 45], these building blocks can be expressed as 



J;'Mp) = T;^{k,p){T^,T'}+T-^{k,p) 
Jt{Kp) = f+^ik,p){T^,T''} + f-^{k,p) 



rj-ia rj-il 



(36) 
(37) 



where the first part is symmetric, while the second one is antisymmetric with respect to the 
color structures of the two gluons. Here k{p), a{b), and fi^u) are the momentum, color, and 
polarization of the off-shell gluons. Below, we refer to the gluon with indices {u, b,p) as the one 
connecting with the spectator quark from the B meson. 
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g'(n,a,k) 




g'(v. b, p) 




JZ{k,p) 



Figure 8: Building blocks Jnt{k,p) (associated with operators Qi^s) and J°;^{k,p) (associated with 



operators Q^fi) for Figs. 6(d) and 6(e). 



In the NDR scheme, after integrating over the (shifted) loop momentum, we can present 
the quantities T^j^{k,p) and T^^^^k.p) as [44, 45] 
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/e ■ p /c ■ p k ■ p 



Aii9 + 



k^Pu 
k ■ p 



T^u{k,p) = aT^y{k,p), 

ik,p) = T-^{k,p) +^ I^g^ 



k ■ p 
4 



k ■ p 



where the matrix E in Eq. (38) is defined by 
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1-75) 



(I-75), (39) 
(40) 

(41) 



(42) 



with the second fine obtained in a four dimension context with the Bjorken-Drell conventions. 
The parameter a in Eq. (40) denotes the chiral structure of the local four-quark operators in 
the weak interaction vertex with a = ± corresponding to {V — A) ^ {V ^ A), respectively. The 
dimensionally regularized expressions for the Ai functions are collected in Appendix B. 

Equipped with the explicit form for these building blocks, we can now evaluate all the 
Feynman diagrams in Fig. 6. After direct calculations, the final results with the subscript 
denoting the contraction of the corresponding operator in the weak interaction vertex are 
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(49) 



8 V(l - n) (1 - v) V ^ 

where the argument is the quark mass propagating in the fermion loops. At this stage, the 
Aj functions are the ones that have been performed the Feynman parameter integrals, whose 
explicit forms can be found in Appendix B. 

With the individual operator contributions given above, the total contributions of these 
higher order h Dg*g* strong penguin diagrams to the decay amplitudes of i? ^ MiMi 
modes can be written as 

Gp 



A, 



b^Dg'g* 



V2 



Cg^ A„ + C\ Aq^ + Ca Aq.^ + C4 Aq^ + Cq Ac^ 



(50) 



In order to specialize these general results to 5 ^ vr Ji , tttt decays, we just need to replace Mi 
and M2 with the corresponding mesons. Explicitly, the h Dg*g* strong penguin contributions 
to the decay amplitudes of the four B irK and the three B ^ tttt decay channels are 



A'{B' 
V2A'{B- 



TT 
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TT 



A'{W ^ TT+A-) = Ab^sg^g* (Ml ^7r,M2^K), 

-V2A'{B^ ^ ir^lf) = Ab^sg*g* (Ml ^ vr, M2 ^ i^) , (51) 
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A'{B~^n-Tr^) = 0, (52) 

where the superscript '/' is indicated there to be distinguished from the next-to- leading order 
results given by Eqs. (18) and (19). The total decay amplitudes are then the sum of these two 
pieces. 

With the total decay amplitudes, the branching ratio for B M1M2 decays reads 

B(B M1M2) = A(B M1M2) +X(B ^ M1M2) ^ ■ S, (53) 

where tb is the lifetime of the B meson, S = 1/2 if Mi and M2 are identical, and 5 = 1 
otherwise, pc is the magnitude of the momentum of the final-state particle Mi^2 in the B meson 
rest frame and given by 

Pc = 7r—\/['^B - ("^A/i + [^B - ("^Mi - rnM^y] ■ (54) 

As for the direct CP asymmetries, we use the definition of the difference of the i?-meson 
minus i?-meson decay rates divided by their sum. With the branching ratios of the CP- 
conjugated modes denoted by B{B /), the CP-averaged branching ratios and the direct CP 
asymmetries for P ^ / decays can be expressed respectively as 

1 



B 

2 



B{B^ f)+B{B^ f) , (55) 



B{B J) - BjB f) 
"^^^ - BiB^f) + BiB^f)- ^^^^ 



3 Numerical calculation and Discussions 
3.1 Input parameters 

The theoretical predictions with the QCDF approach depend on many input parameters such 
as the CKM matrix elements, Wilson coefficients, hadronic parameters, and so on. We present 
all the relevant input parameters as follows. 

-Wilson coefRcients. The Wilson coefficients Ci{fi) in the effective weak Hamiltonian have 
been reliably evaluated to the next-to-leading logarithmic order. To proceed, we use the fol- 
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lowing numerical values at = mb scale, which have been obtained in the NDR scheme [49, 58] 

Ci = 1.082, C2 = -0.185, C3 = 0.014, C4 = -0.035, 

C5 = 0.009, Ce = -0.041, Cj/a = -0.011, Cg/a = 0.059, (57) 

'77 



Cg/a = -1.241, Cio/a = 0.218, C^^ = -0.299, = -0.143. 



-The CKM matrix elements. The widely used parametrization of the CKM matrix ele- 
ments in analyzing i?-meson decays is the Wolfenstein parametrization, which emphasizes the 
hierarchies among its elements and is expanded as a power series in the parameter A = \Vus\ [59], 

( 1-T ^ A\%p-ir])\ 



VcKM 



2 

-A 1 - f 



C(A^) . (58) 



The values of the four Wolfenstein parameters (A, A, p, and 77) could be determined from the 
best knowledge of the experimental and theoretical inputs. In this paper, we take 

A = 0.8533, A = 0.2200, p = 0.20, r/ = 0.33, (59) 

as our default input values [60]. The parameters p and f/ are defined byp = p(l — ^),f/ = 
^ (!-¥)• 

-Masses and lifetimes. For the quark mass, there are two different classes appearing in 
the QCDF approach. One type is the pole quark mass which appears in the evaluation of the 
penguin loop corrections, and is denoted by rriq with q = u,d, s,c,b. In this paper, we take 

mu = md = ms = 0, = 1.47 GeV, = 4.80 GeV, (60) 

as our default input values. 

The other one is the current quark mass which appears in the equations of motion and is 
used to calculate the matrix elements of the penguin operators as well as the chiral enhancement 
factors . This kind of quark mass is scale dependent. To get the corresponding value at 
the given scale, we should use the renormalization group equation to run them, which can be 
found, for example, in [49]. Following Ref. [25], we hold (m„ + md)/fris fixed and use as an 
input parameters. Explicitly, we take 

m„(2 GeV) = md(2 GeV) = 0.0413m,(2 GeV) , 

m,(2GeV) = 90MeV, mb(mfe) = 4.40 GeV . (61) 

21 



where the difference between the u and d quark is not distinguished. 

For meson masses and the hfetimes of the B meson, we adopt the center values given by [60] 



tb^ = 1.671 ps, tb, = 1.536 ps, m^^ = 5.2794 GeV, m^, = 5.2790 GeV, 
mK± = 493.7 MeV , m^fo = 497.6 MeV , m^± = 139.6 MeV , m^o = 135.0 MeV . 
-Light-cone distribution amplitudes of mesons. Since the QCDF approach is based on 
the heavy quark assumption, to a very good approximation, we can use the asymptotic form 
of the LCDAs for light mesons [53, 56, 61] 



With respect to the endpoint divergence associated with the momentum fraction integral over 
the LCDAs appearing in this paper, in analogy to the treatment in Refs. [21, 62], we regulate 
the integral with an ad-hoc cut-off 



with Ah = 500 MeV, and do not distinguish whether this divergence comes from the hard 
spectator rescattering or from the annihilation contributions. The possible complex phase 
associated with this integral has also been neglected. 

As for the B meson wave functions, within our approximation, we need only consider the 
first inverse moment of the LCDA $f (0 defined by [21] 



where the hadronic parameter has been introduced to parameterize this integral. This 
parameter has been evaluated using different methods [63, 64] recently. In this paper, we take 
Xb = 460 MeV as our default input value [63] . 

-Decay constants and transition form factors. The decay constant and the form factors 
are nonperturbative parameters and can be determined from experiments and/or theoretical 
estimations. For the decay constants, we take 



$m(x) =6x(l-x), $p(x) = l. 



(62) 




(63) 




(64) 



fB = 200 MeV [25] 



fK = 160 MeV , 



TT 



131 MeV. 



(65) 



For the form factors involving the B ^ K and B ^ tt transitions, we take 



Fo^^'^(O) = 0.258 







(0) = 0.331 



(66) 
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as the default values at the maximum recoil. In addition, we use the formula 

to parameterize the dependence of the form factor on the momentum-transfer g^, with the fit 
parameters given by 

4 = 0.258, mj.j(7r) = 33.81, rf = 0.330 , m%{K)=37A6. (68) 
All of these values are taken from the latest QCD sum rule analysis [65]. 

3.2 Numerical results and discussions 

With the theoretical expressions and the input parameters given above, we can now evaluate 
the branching ratios and the direct CP asymmetries for B ttK and 5 — tttt decays. For 
each quantity, we first give the predictions at the next-to- leading order in as, and then take 
into account the b Dg*g* strong penguin corrections, which are of order a^. The combining 
contributions of the two pieces, denoted by 0{as+<yl), is then given in the last. For comparison, 
the NF results are also presented. All the averaged experimental data are taken from HFAG [28]. 

3.2.1 The CP-averaged branching ratios for B vr/i , tttt decays 

In the SM, the four B ttK decays are dominated by the b s strong penguin diagrams, with 
additional subdominant contributions from the tree and electro-weak penguin diagrams. The 
three S — vrvr decays, however, are tree-dominated modes. It is therefore expected that these 
higher order strong penguin diagrams considered in this paper should contribute effectively to 
B TxK modes, while have only a minor impact on i? ^ tttt ones. Numerical results of the 
CP-averaged branching ratios for these modes are collected in Table 1. 

The dependence of these CP-averaged branching ratios on the weak phase 7 is shown by 
Fig. 9 (without the annihilation contributions) and Fig. 10 (with the annihilation contributions), 
where the solid and dashed curves correspond to the theoretical predictions with and without 
the b Dg*g* strong penguin contributions included, respectively. The horizontal solid lines 
denote the experimental data as given in Table 1, with the thicker one denoting its center value 
and the thinner ones its error bars. In these and the following figures, the default values of all 
inputs parameters except for the CKM angle 7 are used. 
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Table 1: The CP-averaged branching ratios (in units of 10~^) for B ttK, txti decays with the 
default input parameters. and B^^'^ denote the results without and with the annihilation 
contributions, respectively. The NF results, which are of order are also shown for 

comparison, p = 0.20 and f] = 0.33 



Decay Mode NF Exp. 



B- - 




10.07 


13.28 


17.31 


16.04 


20.44 


24.1 ± 1.3 


B- - 




5.69 


7.30 


9.37 


8.72 


10.97 


12.1 ±0.8 




TX+K^ 


7.71 


10.25 


13.61 


12.46 


16.15 


18.2 ±0.8 






3.38 


4.63 


6.26 


5.70 


7.50 


11.5 ± 1.0 




TT+TT^ 


7.41 


7.69 


7.99 


8.32 


8.63 


4.5 ±0.4 


B- - 




5.12 


5.06 


5.06 






5.5 ±0.6 


b"- 




0.15 


0.16 


0.19 


0.17 


0.21 


1.45 ±0.29 



From these two figures and the numerical results given by Table 1, we can see that: 

• For penguin-dominated B ttK decays, due to the enhancement of the penguin am- 
plitudes, the QCDF scheme prefers larger branching ratios than the NF approximation. 
With our default input parameters, however, predictions for the branching rations are 
still smaller than the experimental data even after the inclusion of the annihilation con- 
tributions, if we consider only contributions up to the next-to-leading order in a^. The 
effects of these higher order b sg*g* strong penguin corrections are very prominent 
in these penguin-dominated B ttK decays. With our input parameters, we find that 
these higher order strong penguin contributions can give ~ 30% enhancement to the 
corresponding branching ratios, and such an enhancement can improve the consistency 
between the theoretical predictions and the experimental data significantly. In addition, 
we find that the effect of the annihilation contributions on the branching ratios, though 
not negligible, is not so large as claimed by pQCD method [2, 51]. 

• For tree-dominated B ^ nir decays, the higher order b dg*g* contributions play only 
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Figure 9: The 7 dependence of the CP-averaged branching ratios for B ttK, txtx decays 
without the annihilation contributions. The sohd and dashed curves correspond to the the- 
oretical predictions with and without the h Dg*g* strong penguin contributions included, 
respectively. The horizontal solid lines denote the experimental data as given in Table 1, with 
the thicker ones being its center values and the thinner its error bars. 



a minor role. To a very good approximation, the B — > vr vr" decay can be considered 
as a pure tree process, and it does not receive annihilation contributions either. The 
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Figure 10: The 7 dependence of the CP- averaged branching ratios for B ttK, -kit decays 
with the annihilation contributions included. The meaning of the curves and the horizontal 
solid lines is the same as in Fig. 9. 



theoretical prediction for the corresponding branching ratio agrees with the data quite 
well. For the other two B ^ -ktx modes, however, theoretical predictions with QCDF 
approach are quite inconsistent with the measured ratios, even with the annihilation 
and the higher order strong penguin contributions included. With our input parameters, 
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we find that the theoretical prediction for B*^ tt^tt^ mode is about an eighth of the 
experimental data; for tt+tt^ mode, on the other hand, a value about two times 

larger than the data is predicted. 

• As for the 7 dependence of the corresponding branching ratios, we can see that the two 
decay modes, B^ tt^tt^ and B^ ir^K'^, are almost independent of this angle, 
since the corresponding decay amplitudes have to a good approximation only a single 
weak phase. In addition, the discrepancy between the theoretical prediction and the 
experimental data for B^ tt+tt" can be removed if we use a large angle 7 ~ 120°. With 
the annihilation and the higher order strong penguin contributions included, the four 
B TiK modes, however, prefer a smaller value for this angle around 7 ~ 80°, which is 
quite consistent with the latest direct experimental measurement 7 = 81° ± 19°(stat.) ± 
13°(si/s.) ± irimodel) [66]. 

• The theoretical predictions for the branching ratios are very sensitive to the value of the 
form factor F^'^'^. For example, the large measured decay rates for the four B -kK 
decays can be well accommodated with a larger value of the form factor as shown by 
Beneke and Neubert [25]. On the other hand, the prediction for B^ t\^t\~ decays can 
become consistent with the data only when a smaller value is used. The large measured 
ratio for 5*^ tt^tt*^, however, remains unresolved with the varying of these parameters. 
It is a tough theoretical challenge to accommodate the current experimental data in the 
SM. 

Since the uncertainties in the predictions for branching ratios can be largely eliminated by 
taking ratios between them, we now discuss the variations of the quantities defined by Eqs. (2)- 
(6) with the higher order h Dg*g* strong penguin contributions included. It is the known 
"vrii'" puzzle [29, 30] that the SM predictions are inconsistent with current experiment data. 
The theoretical predictions and the current experimental data for these ratios are collected in 
Table 2. For the 7 dependence of these quantities, we display them in Fig. 11, where the curves 
and the horizontal solid lines have the same interpretations as in Fig. 9. 

From Table 2 and Fig. 11, we can find that the two ratios Rc and i?„ are indeed approx- 
imately equal within the SM as claimed in Ref. [30], while the experimental data for the two 
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Table 2: Ratios between the CP- averaged branching fractions for B txK, txtx modes. The 
values in the parentheses are the ones without the annihilation contributions. 





NF 


0{a,) 


0{as + al) 


Exp. 




1.272 


1.119 (1.209) 


1.077 (1.163) 


2.20 ±0.31 




0.040 


0.041 (0.042) 


0.048 (0.047) 


0.67 ±0.14 


R 


0.833 


0.845 (0.840) 


0.860 (0.855) 


0.82 ±0.06 


Rc 


1.130 


1.087 (1.100) 


1.074 (1.083) 


1.00 ±0.09 


Rn 


1.140 


1.092 (1.106) 


1.077 (1.087) 


0.79 ±0.08 
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Figure 11: Ratios of the CP-averaged branching fractions defined by Eqs. (2)-(6) as functions 
of the weak phase 7. The meaning of the curves and the horizontal solid lines is the same as 
in Fig. 9. 

quantities are quite different with the puzzling pattern P„ < 1. On the other hand, the value of 
the quantity R predicted by the QCDF approach is well consistent with the experimental data. 
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For the other two ratios and Rqq, the discrepancies between the theoretical predictions 
and the experimental data are quite large. As the b Dg*g* strong penguin contributions 
to B ^ irK, TTTT decays are similar in nature, and hence eliminated in the ratios between the 
corresponding branching fractions, the patterns of the these quantities remain unaffected even 
with these new strong penguin contributions included. From the 7 dependence of the ratios 
between the four B irK decays, a smaller value for this phase is preferred. On the other 
hand, a larger value for this phase is favored hj B ^ titi decays. These inconsistences may 
be hints for new physics playing in the electroweak penguin sector as suggested by Buras et 
al. [30]. 

3.2.2 The direct CP asymmetries for B nK, tttt decays 

Contrary to the NF approximation, the QCDF scheme can predict the strong interaction phases 
and hence the direct CP asymmetries in the heavy quark limit. The numerical results and the 
experimental data for this quantity involving the four irK and the three tttt final states are col- 
lected in Table 3. The 7 dependence of the direct CP asymmetries is displayed in Fig. 12 (with- 
out the annihilation contributions) and Fig. 13 (with the annihilation contributions), in which 
the curves and the horizontal solid lines also have the same interpretation as in Fig. 9. 
From these two figures and the numerical results given in Table 3, we can see that: 

• The direct CP asymmetries for B nK, tttt decays are predicted to be typically small 
with the QCDF formalism. This could be well understood, since the direct CP asym- 
metries are proportional to the sines of the strong interaction phases, which are usually 
suppressed by and/or AQco/f^b within the QCDF formalism. Due to a potentially 
large relative phase between the QCD penguins and the coefficient 02, the fi^ tt^tt" 
mode, however, is an exception to this general rule. The direct CP asymmetries for this 
mode is predicted to be about 55%. 

• Although the individual Feynman diagram in Fig. 6 carries large strong phase, the com- 
bining contributions of these b Dg*g* strong penguin diagrams contain only a relatively 
small one. Thus, these higher order strong penguin contributions to the direct CP asym- 
metries are also small. 
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Table 3: The direct CP asymmetries (in units of 10~^) for B ttK, tttt decays with the 
default input parameters. A(jp and denote the results without and with the annihilation 
contributions, respectively. 



CP CP 

Decay Mode Exp. 



B- - 




0.73 


0.52 


0.65 


0.46 


-2.0 ±3.4 


B- - 




7.59 


6.94 


6.56 


6.07 


4±4 


b'- 




5.31 


4.83 


4.39 


4.08 


-10.9 ± 1.9 


b'- 




-3.08 


-2.84 


-2.71 


-2.54 


-9± 14 


b"- 




-4.73 


-5.51 


-4.54 


-5.27 


37± 10 


B- - 




-0.30 


-0.31 






-2 ±7 


b'- 




55.52 


58.53 


55.03 


55.50 


28 ±39 



• The theoretical predictions for Acp{B vr^vr") and Acp{B 7r~^K~) are quite 
smaller than the experimental data, particularly with the opposite sign. How to accom- 
modate these discrepancies in the SM is still a challenge. 

4 Conclusions 

In this paper, we have revisited the B -kK, tttt decays in the framework of QCDF with the 
b Dg*g* strong penguin contributions included. The main conclusions of this paper are: 

1. For penguin-dominated B txK decays, the higher order strong penguin contributions 
induced by 6 ^ sg*g* transitions to the branching ratios are rather large. With our 
input parameters, we find that these higher order strong penguin contributions can give 
~ 30% enhancement to the corresponding branching ratios, and such an enhancement can 
improve the consistency between the theoretical predictions and the experimental data 
significantly . 

2. For tree-dominated B ^ -k-k decays, the higher order h dg*g* contributions to the 
corresponding branching ratios are quite small. 

30 



25 50 75 100 125 150 175 
y (deg) 



25 50 75 100 125 150 175 
y (deg) 



2 . 



-2 . 



5 
. 5 

. 5 

-5 
-7 . 5 
-10 
-12 . 5 



25 50 



75 100 125 150 175 
y (deg) 



5 


-5 
-10 



™ -15 
^ -20 



25 50 



75 100 125 150 175 
Y(deg) 




25 50 75 100 125 150 175 
y (deg) 



25 50 75 100 125 150 175 
y (deg) 



Figure 12: The 7 dependence of the CP asymmetries without the annihilation contributions. 
The meaning of the curves and the horizontal solid lines is the same as in Fig. 9. 



3. Because of large cancellations among the b — > Dg*g* strong penguin contributions, only 
a relatively small strong phase is remained, So that the contributions have small effects 
on predictions of the direct CP asymmetries. 

4. Since corrections of these higher order strong penguin diagrams to the decay amplitudes 
are similar in nature, and hence cancelled in the ratios between the corresponding branch- 
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Figure 13: The 7 dependence of the CP asymmetries with the annihilation contributions 
included. The meaning of the curves and the horizontal solid lines is the same as in Fig. 9. 



ing fractions, the patterns of the quantities R, Rc, Rn, R^ , and i?oo remain unaffected 

compared to the next-to-leading order results. So we haven't found solution to the "vr-ft'" 
puzzle. Our results indicate that to resolve the puzzle we may have to resort to new 
physics contributions through the electroweak penguin sector as observed by Buras et 
al. [30]. 



32 



5. The theoretical predictions for the branching ratios and tlie direct CP asymmetries still 
have large theoretical uncertainties. The dominant errors are induced by the uncertainties 
of the Fq^""'^ [q^) form factors, strange quark mass ms{iJ,), and the CKM angle 7. 

Although the results presented here have still large theoretical uncertainties, the b Dg*g* 
strong penguin contributions to two-body hadronic 5-meson decays, particularly to penguin- 
dominated modes, have been shown to be very important. Further systematic studies on these 
higher order contributions to charmless B decays are therefore interesting and deserving. 

Note added: After this work is finished, we note an interesting study of corrections to 
B ttK, tttt decays has been carried out by Li, Mishima and Sanda [67] in PQCD formalism. 
However, the contributions studied here as depicted by the Feynman diagrams in Figs. 5 and 6 
are not included in their paper. 
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Appendix A: Correction functions at next-to-leading or- 
der in as 

In this appendix, we present the explicit form for the correction functions appearing in the 
parameters and bi. Details about the calculation can be found in Refs. [21, 25]. 
-One-loop vertex corrections. The vertex parameters Vi{M2) result from the first four 
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diagrams in Fig. 2, given by (with M2 = vr, or K) 



Vi{M2 



du $j\/2 {u 



121n—- 18 + 



12 In 



/i 



6 - fi'(l - m) 



1-4,9,10, 



i = 5,7, 



6,8, 



(69) 



with 



9W 



2u 



1 - u 



In U — ZTT 



(70) 



The scheme-dependent constants —18, 6, —6 are specific to the NDR scheme for 75. ^^^d 

denote the leading- twist and twist-3 LCDAs of the emitted meson M2, respectively. 
-Penguin contractions. The QCD and electro-weak penguin parameters Pfg and -Pg 10 arise 
from the diagrams in Figs. 2(e) and 2(f) . Considering the fact that there exist two distinct 
penguin contractions as shown in Fig. 14, these penguin contributions can be written as 



Pi{M2 



+ C3 



4:, 2 

3ln- + --G,,,(s,) 

mi) 4 



In— + --Gm,(0)-Gm,(1) 

6 jJ, 6 



+ (C4 + Ce) 



1^ In - (n^ - 2) GmM - GmA^c) - Gm,{1] 



1 — u 



+ C3 



4 mfo 2 ^ 
-In — + - - Sp 

3 yU 3 



|ln^ + ^-GM,(0)-GM.(l) 
3 /i 3 



+ (C4 + Ce) 
-2Cf I, 



4n/ rUb 
3 /X 



2) GMM-GM,is,)-G 



(71) 



(72) 



a 



9nN, 



(Ci + N,C2 



4 , m;, 
- m — 

3 yU 



a 



97rNr 



(Ci + iV.Cs) 



4 , m;, 
- m — 
3 n 



Gm^{s^ 



(73) 
(74) 
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Figure 14: Two different penguin contractions. 



where Cp 



2^ , and A'^c = 3 is the number of colors, = 5 is the number of hght quark 
flavors. The pole quark mass ratios, = 0, Sc = (mc/mfe)^, are involved in the evaluation of 
these penguin diagrams. The function Gm2{.s) and Gm2{.s) are defined, respectively, by 

Gm2 {s) = [ du G{s, 1-u) (u) , 
Jo 



Gm^{s) 



duG{s,l-u) 



(75) 



with 



G{s,u) = —A dx X {1 — x)\n[s — x{l — x)u — iS], (76) 
Jo 

where the term i6 is the "e-prescription" . The interpretation of $1/3 ^p^^ is the same as in 
the discussion of vertex corrections. 

-Hard spectator interactions. The parameters Hi{MiM2) originate from the hard gluon ex- 
change between the meson M2 and the spectator quark (corresponding to the last two diagrams 
in Fig. 2) with the results given by 



Hi{M^M2) 



Ib fhli 



«-rnlijK^'''{ml,J Jo ^ 

{l-u){l-v) u{l-v) 



fl rl 

X I du dv 
10 Jo 



(77) 



for i = 1-4,9,10, 



HiiM,M2) 



Jb /mi 



B 



B^Mi , 



m 



M2' 



Jo 4 



»1 A 

X I du I dv 
10 Jo 



u{l-v) {l-u){l-v) 



(78) 



for 2 = 5, 7, and Hi{MiM2) = for i = 6,8. In these results $f (0 is the leading twist LCDAs 
of the B meson as defined by Eq. (14). 
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-Weak annihilation contributions. The basic building blocks for annihilation contributions 
originate from Fig. 3 and given by (omitting the argument M1M2 for brevity) 



A\ 
A{ 

At 



0, 

Txas 

0, 



dudv { $A/2 [u] $ 



+ 



f (1 — uv) u^v 



Ml 



Txas / dudv rf^<l>A/,(n)$;"^(tO 



u{l — uv) 
2v 



+ 



uv^ 



.Ml 



UV 



uv 



dudv < r^^^ ^ivhiu) ^p'[v 



uv{l — uv) ^ 
2(l + u) 



M,{v) ^:-(u 



- p \ 



2u 



■p 



uv{l — uv) 
2(l + v)] 



(79) 



where the superscripts 'z' and '/' refer to gluon emission from the initial and final-state quarks, 
respectively. The subscript '/c' refers to one of the three possible Dirac structures Fi ® F2, i.e., 
= 1 for iy -A)®{y -A), A; = 2 for (V - A) ® (^ + A), and A; = 3 for (-2)(5-P) ® (^ + P). 
Considering the off-shellness of the gluon in Fig. 2 and Fig. 3, it is reasonable to evaluate 
the vertex and penguin corrections at the scale [i ~ m^, while the hard spectator scattering 



and the weak annihilations contributions at the scale \ih = y/AhJi with Ah = 0.5 GeV. 



Appendix B: Analytic expressions for the functions 



In the NDR scheme, after performing the loop-momentum integration, we can present the 
analytic expressions for the Aj functions appearing in Eqs. (38) and (39) as 



AZ5 



Ai23 
AZ24 



-F(-)(4 7r/i 



, 1 rl~x 

/ dx / 

Jo 



dyC-^-^ 2mle-2mlxe-2Px^e + 2Px^e 
2p ye — 2p y e + 2p xy e + A{k ■ p) xy e — A{k ■ p) x ye 



r(^)(4 7r^ 



2 



+4C -12Cx -4Cxe] , 

1 rl — x ^ p 

dx dyC"'^"^ 2mle + 2k^xe-2k'^x^e-2mlye 
Jo 

2„2„., r, „2 „.2 , , o „2 „,3 ^ , /I „N „„. , a „N™„.2 



(80) 



g^, ^ ^ ^ ^ " ■■"q . 

2Jt^ x^ y e — 2p^ y^ e + 2p^ y^ e + A{k ■ p) xy e — A{k ■ p) xy^ e 
+AC -12Cy -ACye] , (81) 
^'^'^^ /'^^ /■'""^..r^-i-^^.. (82) 



4F(-) (4 7r/i2)2 (A;.p)e / dx / dyC-^-^xy 
2 Jo Jo 



4rC-){4nfi')Hk-p)e f dx ('"^ dyC^^~i y [l - y) 
z Jo Jo 



(83) 
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4r(-) (47r/i2)2 (A;.p)e / dx dyC'^'^ x{l- x 



JQ 



1-x 



-AZ23 



2 



r(-)(4 7r^^) 



dx 
Jo 



+2p ye — 2p y e + 2p xy e + A{k ■ p)xy e — A{k ■ p)x ye 
+4C -ACx - ACe + ACxe] , 
r(-) (4 7r/i2)t Tc/x ""dyC-^-^ hmle + 2 x e - 2 k'^ x^ e - 2m^ y e 



JO 

+2 /c^ ?/ e — 2 e + 2 e + 4 (/;; ■ p) X ?/ e — 4 (/c ■ p) X 

+AC -ACy -ACe + ACye] , 

-r(-) (4 7r/i2)f /^c/x ''dyC-^-^ [2 e + 2 x e + 2 fc^ x^ e - 2 P x^ e 
2 Jo Jo I >j <j 

+2p^ y e — 2p^ y"^ e — 2p^ xy'^ e + A{k ■ p) x^ y e 
+4C + ACx-ACe-ACxe] , 
-r(-) (4 7r/i2)t /^c/x ^dyC~^-^ [2 e + 2 P x e - 2 /c^ x^ e - 2 y e 

2 JO JO L « « 

-4/c^x?/e + 2 A;^x^ ye - 2p^y^e + 2p^y^e 

—A{k-p)ye + A{k-p)xye + A{k-p)y'^e — A{k-p)xy'^e 

+AC -ACy-ACe + ACye] , 

h e - 2 X e - 2 x^ e + 2 A;^ x^ e 
2 p'^ y e — Ap^ X y e — 2 p^ y'^ e + 2 p^ X y"^ e 
A{k ■ p) X e + A{k ■ p) x^ e + A{k ■ p) xy e — A{k ■ p) x^ y e 



r(-)(4vr/.2)^ 



dx 
Jo 
2 



+4C-4CX -4Ce + 4Cxe] , 

r(-) (47r/i2)f /^rfx ""rfyC-i-t [2m^e + 2A:2xe-2A:2x2e + 2m2?/e 
2 Jo Jo L « « 

-2 /c^ x^ ?/ e + 2 e - 2 e + 4 (A; ■ p) X 2/^ e 
+4C + 4C|/-4Ce-4C|/e] , 

l-x 



r(-) (4 7r/i2)5 (A; ■p)e / c/x / c/yC-i-f (1 - x) x 



JO 

1 



4r(-) (4 7r/i2)2 (A; .p)e / dx dy C'^-^ x {I - x) {I - 2y) 



JO 

1 



l-x 



-AT{-){ATi^i')i {k-p)e I dx I dyC~^^i xy{l -2x) 



JO 

1 



Jo 



l-x 



AT{-){A'K^^)^ {k-p)e dx dyC~'^~2y{l-x-y + 2xy) 
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Aii9 = 4.V{-){Atiii'^)^ {k-p)e j dx I dyC'^'^ x {1 - x - y + 2xy) , (96) 
2 JO Jo 

A^2o = -4r{U{47rfi')Hk-p)e f'dx f'~'dyC~'-Hl-2x)y{l-y), (97) 
2 JO JO 

A^2i = 4r{^){A7rfi')Hk-p)e f'dx f'"^dyC-'~^^xyil-2y), (98) 

2 JO JO 

A222 = -8r(^)(4V)f (^.p) r^^^ f^-^^^^-l-f (1_^)^2^ (gg) 

2 JO JO 

where the parameter C is defined by 

C = ml - X {1 - x) k'^ - y {1 - y) p'^ - 2 xy {k ■ p) - i 6. (100) 

with mg being the quark mass in the Fermion loops. 

For B meson decaying into two hght energetic hadronic final states, the characteristic scale 
for the quark momentum of the final-state meson constituents is of order m^, whereas the 
momentum of the spectator quark from the B meson is of order Aqcd- Assuming that the 
off-shell gluon with index {u,b,p) is connected with the spectator quark in the B meson, at 
leading power in AQCD/f^ib, the Ai functions given above can then be simplified greatly. After 
subtracting the regulator e using the MS scheme and performing the Feynman parameter inte- 
grals, we get (here we give only the relevant Ai functions needed in this paper; details for the 
others can be found in Ref. [47]) 

At, = 2+'^[Go{n)-Go{n + rs)]--[G^i{n)-G^i{n + rs)], (101) 

A- o 4 2ri(l+r3) 2 (ri + rs + ri ra) 

Ate = -2 \ 2 Go{ri} k Go{ri + rs) 

+ -[G„i(n)-G„i(n + r3)] - ^^'^^ Un) 

^(4-n-r3Hn±lj)j,^(,. (102) 

^^3 

Az23 = -2 - ^ [Go(ri) - Goin + rg) ] + - [G^i{n) - G^,{n + r,)], (103) 

^3 ^^3 
A226 = -A223, (104) 

9 6 rric 3r3 3r3 

+^[G_i(ri)-G_i(ri + r3)] , (105) 



22 12 4ri 8 /i 2 (7ri - r3 - 3ri rg + 2r^ - 2rl; ^ ^ 
9 r3 3r3 3 rric of 3 
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2ri(7 + 2ri-3r3) ^ , 4j2r^+r3) r ^ ^ i 

-J '3 '3 

3(4-ri)ri 3 (4 - n - rg) (ri + rg) 
+ 2 ^o(n) 2 T'ol'-i +^^3) , (106) 

A- 32 16 u 8(2 + ri) ^ , 8 (2 + n + rg) ^ , 



32 12 4 n 16 /i , 2ri(7 + 2ri + 6r3) ^ , , 
9 rg 3r3 3 3r3 
2 (2r2-r3 (l-4r3)+ri (7 + 6r3)) 



3ri 



Go(ri +r3) 



[G_i(ri) - G_i(ri + r3) ] + ^ 2'"'^''' ^o(ri) 

^3 n 

3 (4 - ri - rs) (ri + rs) 



To(ri+r3) (108) 
'3 

2 2 (8 + ri) ^ , , 2 /8 + ri 4 \ , , 
A... = 3 + ^37r^^°(^^)-3 (^ + ;7T7;)^°(^^ + ^^) 

-- [G_i(ri)-G_i(ri+r3)] , (109) 

^^3 

2 16 8ri 2ri(4r2 + 3r3(8 + r3)+ri(20 + 7r3)) „ . , ^ 
A221 = — + o 2 / , N G'o(ri+r3) 

3 r3 3r3 3r| (n + r3) 

2ri (20 + 4ri + 3r3) ^ , . ^ 4(4ri+r3) , . n I ^ \-\ 

ir2 G'o(n) + 2 G'-i(^i) - <^-i(n + ^3) 

3r| r| 

" r,(n) + ^ <^ - - 7' + ■'^l ro(n + . (no) 

n n 

where we have introduced the notations = k'^/m'^, r2 = p^/m^ and = 2{k ■p)/m'^^, 

with niq = rric or rrib. For hght u, d, s quark loops, these Ai functions can be evaluated 
straightforwardly. 

The functions Gi(t) {i = —1,0) are defined by 

Gi{t)= [^dxx'\n[l-x{l-x)t-i6] . (Ill) 
Jo 

The explicit form for G^ifi{t) could be found in Ref. [45]. 

In addition, we have also introduced the function Tj(t), which is defined by 

^■<'>-/o'''" l-.(l-'.)t-a - '"2' 



The explicit form for TQ{t) is given by [47] 

4 arctan - 



To(t) = { ^ ' _ (113) 



< t < 4 
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